jhifl  report  may 


released  to  0IS« 


RESEARCH  LABORATORY 


•5 

i 


ESTIMATING  METAL  PARTICLE 
COMBUSTION  KINETICS 


BY 


R.  W.  BARTLETT 
J.  N.  ONG,  JR. 

W.  M.  FASSELL.  JR. 
C.  A.  PAPP 


1 


2  JULY  1962 


AERONUTRONIC 

AERONUTRONIC  DIVISION  •  FORD  MOTOR  COMPANY 


FORD  ROAD/NEWPORT  BEACH,  CALIFORNIA 


4 


> 


CLEARINGHOUSE  FOR  FEDERAL  SCIENTIFIC  AND  TECHNICAL  INFORMATION 
Loan  Document  Processing  Record 


Accession  No.  /6/a/ <_T  £>  3 


19^ 


Date  received, 
by  Input 
Date  due  out., 
from  Input 


THIS  DOCUMENT  IS  ON  LOAN 


Please  fill  in  date  and  sign  when  each  action  below  has  been  completed. 


Film  for  *ttcrofiche 


Rate 


Film  for  Prestock  or  Restock 


Report  Returned  to  Input 


Document  Returned  to  Sender 


I 


TS-J03-2/65 


To  CASE  FILE  after  Document  is 
Dispatched 


CLEARINGHOUSE  FOR  FEDERAL  SCIENTIFIC  AND  TECHNICAL  INFORMATION,  CFSTI 

INPUT  SECTION  4IS.II 


LIMITATIONS  IN  REPROOUOTION  QUALITY  OF  TECHNICAL  ARSTRACT  BULLETIN 
DOCUMENTS,  DEFENSE  DOCUMENTATION  CENTER  (DDC) 


□  •. 


M 


AVAILABLE  ONLY  FOR  REFERENCE  USE  AT  DDC  FIELD  SERVICES. 
COPY  IS  NOT  AVAILABLE  FOR  PUBLIC  SALE. 

AVAILABLE  COPY  WILL  NOT  PERMIT  FULLY  LEGIBLE  REPRODUCTION. 
REPRODUCTION  MILL  BE  MADE  IF  REQUESTED  BY  USERS  OF  DDC. 


fyj  A.  COPY  IS  AVAILABLE  FOR  PUBLIC  SALE. 

□  B.  COPY  IS  NOT  AVAILABLE  FOR  PUBLIC  SALE. 

n  3  LIMITED  NUMBER  OF  COPIES  CONTAINING  COLOR  OTHER  THAN  BLACK 

AND  WHITE  ARE  AVAILABLE  UNTIL  STOCK  IS  EXHAUSTED.  REPRODUCTIONS 
WILL  BE  MADE  IN  BLACK  AND  WHITE  ONLY. 


DATE  PROCESSED:  <5“//* /<* ,5- 
PROCESSOR: 


TSL-I2I-2/S4 


.  .  I  wnmifti  WBOt  JgggTJi^  (fWdH*TU'i 


tamtrx  ffWFR»nw9»watiw« 


awawift  iMmwmmmimtmimxw'm 


PUBLICATION  NO.  U-I75V 


RESEARCH  LABORATORY 


ESTIMATING  METAL  PARTICLE 
COMBUSTION  KINETICS 


BY 


R.  W.  BARTLETT 
J.  N.  ONG.  JR. 

W.  M.  FASSELL,  JR. 
C.  A.  PAPP 


2  JULY  1962 


THE  WORK  DESCRIBED  HEREIN  WAS  SPONSORED  BY  THE 
ADVANCED  RESEARCH  PROJECTS  AGENCY  THROUGH  THE 
BUREAU  OF  NAVAL  WEAPONS  UNDER  CONTRACT  NOw  61-0905-e 


AERONUTRONIC  DIVISION  •  FORD  MOTOR  COMPANY 


FORD  ROAD /NEWPORT  BEACH,  CALIFORNIA 


ABSTRACT 


A  combustion  model  for  metal  particles  is  derived  and  compared 
with  experimental  results  for  a  few  metals  with  high  combustion  enthalpies. 
It  is  concluded  that  aluminum  and  beryllium  oxidize  slowly  because  of  for¬ 
mation  of  solid  oxide  films  which  impede  the  transport  of  reactants  at 
temperatures  below  their  oxide  melting  points.  Diffusion  through  the 
liquid  oxides  is  several  decades  faster,  causing  ignition  to  occur  at  the 
oxide  melting  point.  After  ignition,  the  temperature  rises  rapidly.  When 
the  metal  boiling  point  is  reached,  metal  evaporation  expands  the  liquid 
oxide  shell  until  the  evaporation  rate  and  rate  of  diffusion  of  reactant 
through  the  oxide  shell  are  equal.  The  oxidation  reaction  occurs  at  the 
oxide  shell  rather  than  the  metal  interface  and  proceeds  by  this  mechanism 
until  the  particle  is  consumed.  Expansion  of  this  shell  accounts  for  the 
prevalence  of  hollow  oxide  spheres  observed  in  the  quenched  combustion 
products.  Calculations  of  the  elapsed  times  for  each  stage  in  the  com¬ 
bustion  sequence  have  been  machine  computed .  Inadequate  knowledge  about 
diffusion  coefficients  in  liquid  oxides  is  the  principal  uncertainty  in¬ 
volved.  The  diffusion  flame  theory  adequately  accounts  for  the  combustion 
of  magnesium. 
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C  ■  reactant  concentration  (g  cm  ) 

AC  »  concentration  difference  across  diffusion 
barrier  (g  cm“3) 

C  *  heat  capacity  at  constant  pressure,  per 
**  unit  mass  (cal/g°K) 


effective  gas  molecule  collision  diameter 
(cm) 


Dq  ■  initial  particle  diameter  (cm) 

D12  -  diffusion  coefficient  of  reactant  species 

in  liquid  oxide  (cm^  sec"^-) 

AF^  =  activation  free  energy  (kcal/mole) 

g  *  acceleration  of  gravity  (980  cm  sec“^) 

h  =  heat  transfer  coefficient  (cal  cm"^  sec"^  °K"^) 


h  =  Planck's  constant 


k  =  Boltzmann  constant 

k  *  thermal  conductivity  (cal  cm-^  s~V  °K“^) 

k  »  conversion  constant 

o 

kl  =  Wr  ®p  ^  nkv 

K  =  ratio  of  R  /R. 

m  i 

m  *  mass  of  metal  oxidized  (g) 

M  ■  molecular  weight 

q  ■  heat  flux  (cal  cm-4  sec"A) 
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total  rate  of  heat  flow  to  metal  sphere 
from  metal  vapor  (cal  sec"^) 

radius  of  isothermal  contours  (cm) 

gas  constant 

particle  radius  (cm) 

inner  radius  of  oxide  shell  (cm) 

outer  radius  of  oxide  shell  (cm) 

time  (sec) 

temperature  (°K) 

hot  gas  streaming  velocity  (cm  sec-^) 
relative  to  particle 

metal  vapor  diffusion  velocity,  radial 
direction 

mass  transfer  rate  — f  (g  sec  ) 

at 

metal  boiling  point  (°K) 

2  -1 

thermal  diffusivity  (cm  s  ) 

latent  heat  of  metal  fusion  (cal/g) 

latent  heat  of  metal  vaporization  (cal/g) 

enthalpy  of  metal  at  ignition  temperature 
(cal/g) 

enthalpy  of  liquid  metal  at  boiling  point 
(cal/g) 

heat  of  formation  of  metal  oxide  (cal  g  ^ 
of  metal) 

gas  viscosity  (poise) 
metal  density  (g  cm”'*) 
oxide  density 
Stefan  Boltzmann  constant 
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Subscripts 


c  =  metal  oxide 

f  =  gas  film  adjacent  to  metal  particle 

g  =  gas  or  flame 

i  ■  ignition 

m  *  metal  or  metal  surface 

o  ■  initial  condition 

v  ■  metal  vapor 

v  ■  radiation 
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INTRODUCTION 


Most  of  the  studies  of  heterogeneous  combustion  have  assumed 
that  the  combustion  reaction  occurs  in  the  gas  phase,  forming  a  spherical 

diffusion  flame  around  the  fuel  droplet.  This  theory  has  been  summarized 

1  2 
by  Penner  and  adequately  accounts  for  the  combustion  of  magnesium. 

However,  experience  with  higher  boiling  metals  has  indicated  that  their 

2-4 

behavior  does  not  conform  to  the  diffusion  flame  theory.  It  has  long 
been  recognized  by  Vulis  and  other  Russian  investigators  that  the  combus¬ 
tion  of  carbon  cannot  occur  by  a  gas  phase  reaction  but  must  occur  at  the 
particle  or  pore  surface.^  They  have  observed  that,  depending  on  the  flow 
and  thermal  regimes,  the  combustion  rate  will  be  controlled  by  either  dif¬ 
fusion  of  oxygen  or  product  gas  species  to  or  away  from  the  surface  or  by 
the  rate  of  chemical  reaction  at  the  surface. 

Only  metals  with  high  oxide  formation  enthalpies  are  attractive 
energy  sources  and  as  a  result  their  oxides  are  invariably  stable  con¬ 
densed  phases  at  high  temperatures.  When  this  situation  is  coupled  with 
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a  moderately  high  metal  boiling  point,  the  slow  and  rate  controlling  step 
may  shift  to  transport  of  reactants  through  the  oxide  layer.  Aluminum  and 
beryllium  are  believed  to  be  examples  of  this  mode  of  combustion.  The 
correlation  of  minimum  aluminum  ignition  temperature  with  the  melting 
point  of  Al^O^  has  been  observed  by  Friedman  and  Macek.^  This  study  is 
principally  concerned  with  the  analysis  of  this  latter  type  of  combustion. 
A  combustion  model  is  presented  and  the  results  calculated  from  it  are 
compared  with  the  available  experimental  data  on  aluminum  and  beryllium. 
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EXPERIMENTAL  SYNOPSIS 


Burning  times  were  determined  by  phc  .ographing  the  incandescent 
particles  in  a  me thane -oxygen  torch  usi-  g  a  high  speed  motion  picture 
camera  or  a  still  plate  camera  in  combination  with  a  rotating  light 
chopper. ^  The  particles  travelled  in  the  same  direction  as  the  flame 
and  ignition  delay  times  were  approximated  by  dividing  the  measured  dis¬ 
tance  of  the  hot  zone  of  the  flame  through  which  the  particle  travelled 
before  igniting,  by  the  particle  velocity.  The  pre-ignition  particle 
velocity  was  assumed  to  be  the  same  as  the  post-ignition  velocity  deter¬ 
mined  by  the  motion  picture  camera  or  chopper.  The  sharp  transition  of 
the  particle  from  nonluminous  to  brilliant  incandescence  was  used  to 
indicate  ignition.  Observations  were  made  on  aluminum,  magnesium, 
beryllium,  zirconium,  titanium,  niobium,  tantalum,  silicon,  chromium, 
iron  and  various  alloys  of  these  metals.  The  principal  qualitative 
results  are  summarized  below.  The  quantitative  results  are  compared 
with  predicted  results  later. 
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(1)  All  of  the  metals  but  beryllium  and  large 
aluminum  and  silicon  particles  ignited  and 
burned  in  the  torch. 

(2)  With  the  exception  of  magnesium,  the  com¬ 
bustion  products  included  as  a  major  con¬ 
stituent  hollow  spheres  about  the  same  size 
as  the  original  particle,  Figure  1.  Pro¬ 
pellent  strand  burning  tests  also  produced 
hollow  spheres. 

(3)  The  MgO  particles  were  much  smaller,  usu¬ 
ally  less  than  one  micron,  and  often  ap¬ 
peared  as  single  crystals.  Magnesium- 
aluminum  alloys  with  a  magnesium  content 
greater  than  20%  also  burn  in  this  manner. 
The  difference  between  magnesium  and  al¬ 
uminum  particle  combustion  is  illustrated 
by  Figure  2,  which  shows  a  mixture  of 
aluminum  and  magnesium  particles  burning 
in  a  torch  flame.  The  luminous  cloud  is 
characteristic  of  magnesium.  The  sharp 
light  spots  are  individual  aluminum 
particles . 
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FIGURE  1 .  CROSS  SECTION  OF  ALUMINUM  ALLOY  COMBUSTION  PRODUCTS 
SHOWING  SPHERICAL  OXIDE  SHELLS 


FIGURE  2.  MIXTURE  OF  -115  +150  MESH  ALUMINUM  AND  MAGNESIUM  PARTICLES 

BURNING  IN  AN  OXYGEN- NATURAL  GAS  FLAME;  PHOTOGRAPHED  THROUGH 
A  160- INTERRUPTIONS-PER- SECOND  LIGHT  CHOPPER 
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(4)  Starring,  attributed  to  fragmentation,  and 
spiral  motion  of  some  metals  and  alloys 


were  often  observed. 

Any  successful  analysis  of  metal  particle  combustion  kinetics 
must  predict  the  burning  times  with  some  reasonable  degree  of  accuracy, 
provide  for  reasonable  agreement  with  observed  ignition  temperatures 
and  explain  the  formation  of  hollow  oxide  spheres. 
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COMBUSTION  MODEL 


Heat  and  mass  transfer  problems  usually  can  be  solved  only  for 
simple  geometric  shapes  and  all  of  the  calculations  presented  here  are 
based  on  spherical  geometry.  The  melting  points  of  most  metals  with  high 
oxidation  enthalpies  are  on  the  lower  end  of  the  combustion  temperature 
range,  300°K  to  >  3000°K.  Since  surface  tension  will  cause  the  liquid 
metal  to  assume  a  spherical  shape,  an  error  due  to  the  assumption  of 
spherical  geometry  for  nonspherical  particles  will  only  exist  during  a 
brief  initial  period. 

Before  presenting  the  method  of  calculation,  a  review  of  the 
model  on  which  it  is  based  is  necessary.  For  a  metal  sphere  at  an  in¬ 
itial  low  temperature  Tq,  injected  into  a  hot  gas  stream  or  flame  at 

temperature  T  ,  the  oxidation  rate  will  depend  on  the  surface  tempera- 
8 

ture  and  mechanism  of  oxidation.  For  many  metals,  including  aluminum 

and  b^rylliu  .,  the  oxidation  rates  are  believed  to  be  controlled  by 

8  9 

diffusion  of  reacting  ions  through  the  solid  oxide  coating  ’  that  is 


formed.  The  rate  is  relatively  slow  as  long  as  the  oxide  remains  solid. 
Diffusion  rates  in  liquids  are  much  faster  than  in  solids.  The  increase 
is  several  orders  of  magnitude^  and  the  oxidation  rate  will  increase 
correspondingly  at  the  oxide  melting  point.  Consequently,  for  this  type 
of  system,  ignition  is  expected  to  occur  at  the  melting  point  of  the 
oxide,  T^. 

The  oxide  coating  formed  on  magnesium  is  not  adherent  and  does 
not  act  as  an  effective  diffusion  barrier.  Oxygen  gas  can  pass  directly 
to  the  metal  surface.  Provided  oxygen  diffusion  through  the  gas  boundary 
layer  is  not  limiting,  the  rate  of  oxidation  of  this  type  of  metal  de¬ 
pends  on  the  surface  reaction  rate  and  increases  exponentially  with 
temperature, 


dm 

dt 


C .  k 
A  o 


-  AF+/RT 


(1) 


As  the  particle  is  heated,  a  temperature  is  eventually  reached  where  the 
thermal  flux  from  oxidation  becomes  appreciable.  This  increases  the 
temperature  at  a  faster  rate  and  very  soon  the  temperature  itself  is 
increasing  at  a  quasi  exponential  rate.  For  this  system  there  will  not 
be  a  precise  ignition  temperature  but  a  narrow  temperature  band  above 
which  the  oxidation  rate  and  temperature  increase  sharply.  An  average 
temperature  in  this  band  is  taken  as  the  ignition  temperature  for  this 
type  of  system.  For  magnesium,  this  has  been  observed  to  be  about  900°K. ^ 
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In  both  cases,  the  heat  flux  and  surface  regression  rate  caused 
by  oxidation  are  negligible  before  ignition  but  must  be  considered  at 
higher  temperatures . 

Following  ignition,  the  combustion  behavior  of  magnesium, 
beryllium  and  aluminum  differ  from  each  other  and  also  depend  on  the 
system  pressure.  Magnesium  has  a  vapor  pressure  of  10  mm  at  1000 °K 
and  its  normal  boiling  point  is  1381°K.  Magnesium  oxide  is  a  stable 
solid  at  these  temperatures  and  since  no  phase  changes  are  involved, 
it  is  expected  to  remain  porous  and  permit  rapid  transport  of  oxygen 
and  also  metal  vapor.  Magnesium  combustion  is  directly  related  ta  the 
vaporization  of  magnesium.  Although  a  direct  gas  phase  reaction  to 
produce  MgO  vapor  is  not  thermodynamically  possible,  magnesium  vapor 
species  and  oxygen  can  react  on  surfaces  to  produce  solid  MgO.  As  a 
result,  very  small  MgO  crystals  quickly  nucleate  and  grow  in  the  sur¬ 
rounding  gas,  accounting  for  the  luminous  cloud  and  small  crystals  that 
are  observed  in  burning  experiments.  The  combustion  rate  is  dependent 
on  the  evaporation  rate  which  depends  on  the  heat  flux  to  the  metal 
particle . 

The  behavior  of  aluminum,  beryllium  and  many  other  metals  is 
quite  different  from  magnesium  even  above  their  ignition  points.  As  a 
result,  longer  periods  are  required  to  complete  combustion  and  an  analy¬ 
sis  of  the  kinetic  process  is  warranted.  Above  the  ignition  temperature, 
aluminum  and  beryllium  are  covered  with  molten  oxide  layers  which  prevent 
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metal  evaporation  provided  the  vapor  pressure  of  the  metal  is  less  than 
the  system  pressure.  The  normal  metal  boiling  point  and  oxide  melting 
point  of  beryllium  are  almost  the  same.  The  temperature  of  aluminum 
must  be  raised  from  2300°K  to  2740°K  for  boiling  to  occur.  In  both 
cases,  once  the  boiling  point  is  attained,  metal  evaporation  will  expand 
the  oxide  shell  providing  a  metal  vapor  envelope  between  the  metal  sphere 
and  the  oxide  shell.  It  is  important  to  emphasize  that  diffusion  through 
the  liquid  oxide  shell  is  still  slow  when  compared  with  gaseous  diffusion 
and  reaction  will  occur  at  the  oxide  shell  rather  than  at  individual 
nuclei  in  the  surrounding  gas.  The  oxide  shell  will  expand  or  contract 
to  permit  the  mass  rate  of  metal  evaporated  from  the  metal  sphere  to 
equal  the  rate  of  oxidation  of  metal  at  the  oxide.  Any  sudden  ejection 
of  the  particle  from  the  hot  zone  quickly  freezes  the  oxide  shell  and 
accounts  for  the  hcvLlow  oxide  spheres  with  contained  metal  particles 
which  have  been  collected,  sectioned  and  observed  microscopically. 

Such  particles  are  very  numerous  in  most  combustion  experiments. 

In  any  practical  thermodynamic  application,  the  chamber  pres¬ 
sure  must  be  much  higher  than  one  atmosphere,  and  as  a  result,  the  metal 
boiling  points  will  increase.  The  extrapolation  of  vapor  pressure  data 
taken  at  temperatures  below  1700°K  to  temperatures  above  the  normal 
boiling  points  is  very  hazardous.  In  lieu  of  better  information,  however, 
such  an  extrapolation  indicates  a  boiling  point  of  about  3450°K  for  both 
aluminum  and  beryllium  at  1000  psi.  As  the  chamber  pressure  and  boiling 
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point  increase,  the  thermal  gradient  causing  evaporation  will  decrease 
and  the  combustion  time  will  increase.  At  sufficiently  high  pressures, 
no  evaporation  can  occur. 

It  is  convenient  to  separate  the  combustion  of  metals  which 
behave  in  the  manner  of  aluminum  and  beryllium  into  three  consecutive 
stages  separated  by  the  ignition  temperature  and  the  metal  boiling 
point  and  compute  the  elapsed  time  for  each  stage.  The  three  stages 
are  shown  schematically  in  Figure  3  and  summarized  in  Table  I.  The  sum 
of  these  elapsed  times  gives  the  time  interval  required  for  complete 
combustion  to  occur.  All  calculations  are  based  on  assumed  flame  tem¬ 
peratures  of  3300°K  for  aluminum,  3450°K  for  beryllium  and  3150°K  for 
magnesium. 
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FIGURE  3.  SCHEMATIC  PLOT  OF  METAL  AND  OXIDE  TEMPERATURES  AND 
METAL  DROPLET  RADIUS  DURING  THE  THREE  COMBUSTION 
STAGES 
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(3)  Oxidation  Diffusion  Co 
trolled  Through  Shell. 


IGNITION  DELAY 


Thermal  Flux  Calculation 

The  thermal  flux  to  the  metal  during  the  period  prior  to 
ignition  is  the  sum  of  convection  and  radiation  from  the  surroundings. 
During  this  period  the  oxidation  reaction  is  a  negligible  thermal  source. 
The  convective  flux,  q^,  across  the  gas  boundary  film  can  be  calculated, 


hf  (Tg  ' 


V 

m 


using  the  Nusselt  number 

Nu 


D 

o 

f 


(2) 

(3) 


The  following  correlations  for  forced  and  free  convection  were  taken 

13 

from  Ranz  and  Marshall.  For  high  streaming  velocities,  the  Nusselt 
number  depends  on  the  Reynolds  and  Prandtl  numbers: 


D 

o 

f 


2.0  +  0.60 


D  v  J , 
o  to  ^f 


1/2 


(4) 
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For  free  convection  the  Nusselt  number  depends  on  the  Grashof  number, 


The  streaming  velocity  of  the  hot  gases  over  a  metal  particle  is  unknown 

i 

and  probably  varies  from  a  low  value  corresponding  to  free  cpnvection  to 
the  flame  velocity  relative  to  the  combustion  chamber.  The  latter  veloc¬ 
ity  is  also  a  variable  quantity  depending  on  the  relative  location  in  the 
flame  and  whether  the  flame  configuration  changes  with  time.  The  thermal 
flux  will  be  lowest  and  the  ignition  delay  longest  for  free  convection 
and  calculations  based  on  the  free  convection  assumption  will  be  more 
conservative  than  arbitrarily  assuming  an  average  streaming  velocity. 

The  diameters  of  metal  particle  additives  are  of  the  order  of  100  microns 

-3 

which  corresponds  to  low  Grashof  numbers,  A1  *  3.37  x  10  ,  Be  *  3.17  x 

-3  -3 

10  and  Mg  =  4.52  x  10  .  As  a  result,  the  second  term  in  equation  (5) 

only  amounts  to  about  10%  of  the  total  and  can  be  neglected.  Therefore, 
for  these  systems,  free  convection  is  approximately  equivalent  to  heat 
transfer  from  a  stagnant  fluid, 


h.D 
f  o 


2 
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and  from  equation  (2) 


qf 


2^ 

D 


(6) 


As  a  consequence,  the  convective  flux  is  inversely  proportional  to  the 
particle  size. 

A  solution  of  equation  (6)  requires  an  evaluation  of  k^.  The 
prediction  of  thermal  conductivities  for  polyatomic  gases  is  only  semi- 
empirical  even  for  low  pressure  systems.  Eucken's  equation  for  thermal 
conductivity  of  polyatomic  gases  is:i 


k 


(7) 


A  good  estimate  of  k  can  be  made  by  employing  the  following  equation 
for  the  viscosity 


H 


2.66  x  10 


(8) 


in  which  d  is  the  collision  diameter  and  0  is  a  dimensionless  function 

varying  slightly  from  unity  and  measures  the  deviation  of  the  gas  from 

rigid  sphere  behavior.  The  values  are  based  on  a  Lennard-Jones  function 

14 

for  the  potential  energy  of  interaction. 


-17- 


For  purposes  of  calculations,  the  primary  combustion  gases  were 
assumed  to  be  CO,  H^O,  and  HC1  in  equal  mole  fractions.  An  average 
film  thermal  conductivity  coefficient  k°  was  calculated  by  determining  the 
thermal  conductivity  coefficient  for  each  of  these  gas  species  using 
equations  (7),  (8)  and: 


■f  -  k!  (I  ) 
i  rl 


i 


Dissociation  of  these  gases  may  be  important  at  high  temperatures  but 
dissociation  was  ignored  for  these  initial  calculations.  The  gas  film 


temperature  was  taken  as: 


T  1/2  +  T  1/2 
R _ m _ 


The  values  of  k°  and  the  initial  and  final  convective  thermal  fluxes 
for  100  micron  spheres  of  magnesium,  aluminum  and  beryllium  are  listed 
in  Table  II. 

The  heat  flux  transferred  by  radiation  to  a  small  sphere  from 
a  large  cavity  is  provided  by  the  Stephen  Boltzmann  law: 

a  -  a  (a  T  4  -  e  T  4)  (11 

mgmm 


The  term  a  is  the  absorptivity  of  the  metal  for  radiation  emanating 
8 

from  T  .  The  initial  and  final  radiant  fluxes  are  also  listed  in  Table 
8 

II.  The  values  are  based  on  estimates  for  effective  emissivity  of  .20 
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TABLE  II 


PREIGNITION  DELAY  DATA  FOR  SPHERES,  D0  =  100  MICRONS 
IN  A  METALLIZED  FLAME  AT  THE  TEMPERATURES  INDICATED 


A1 

Be 

Mg 

=  0.20 

=  0.20 

=  0.20 

1  o  i  -l  Ov“3/2» 

kj  (cal  cm  sec  K  ) 

2.36  x  10-4 

2.44  x  10"4 

2.30  x  10 

qf  (cal  cm  ^  sec 

1  T 

o 

5293 

5870 

4773 

qf  (cal  cm  ^  sec 

Ti 

2432 

1509 

4354 

q  (cal  cm  sec  ) 

T 

o 

32.1 

38.3 

27.0 

,  i  "2  -In 

q  (cal  cm  sec  ) 

Ti 

24.2 

19.2 

25.8 

Elapsed  Time  to  heat  sphere 

from  Tq  to  T^  (milliseconds) 

0.5 

1.32 

0.1 

Elapsed  time  to  fuse  sphere 

(milliseconds) 

0.173 

0.50 

0.12 

Total  preignition  delay 

(milliseconds) 

0.673 

1.82 

0.22 

Estimated  flame  temperature 

(°K) 

3300 

3450 

3150 
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which  is  somewhat  less  than  the  average  values  o£  the  total  emissivities 
for  the  oxides  over  the  preignition  temperature  range. ^  The  calculations 
also  assume  e^  =  a^.  The  estimation  of  high  temperature  emissivities  is 
hazardous,  but  the  radiation  flux  is  less  than  the  convective  flux  and  a 
30%  error  in  e^  will  only  result  in  a  4%  to  6%  error  in  the  final  elapsed 
time  calculation. 

Heat  Conduction  and  Temperature  Increase  Within  The  Particle 

To  determine  the  temperature  increase  in  the  sphere,  the  equation 
of  heat  conduction  in  spheres  must  be  solved.  Because  the  flux  at  the 
surface  is  geometrically  uniform,  isothermal  surfaces  within  the  metal  are 
concentric  spheres  and  the  equation  of  heat  conduction  is: 


ST 

ST 


0!  m 


5 

Sr 


r2  ST 
Sr 


(12) 


The  flux  at  the  surface  is  the  sum:  q  +  q 

r  n  v 


q*  +  q 
f  v 


-k 


m 


ST 

Sr 


r  *  R 


(13) 


The  solution  to  this  problem  is  provided  by  Carslaw  and  Jaeger^  in  the 

kmT  3aml  r  amfc 

form  of  plots  of  — —  -  — r—  versus  —  for  various  values  of  — tt 

Rq  r  R  rz 


for  a  constant  flux.  Calculations  of  the  temperature  lag  at  the  center 
of  100  micron  spheres  when  the  surface  temperature  Is  at  the  Ignition 
point  were  made  assuming  the  maximum  thermal  flux  as  the  constant  flux. 
The  temperature  lag  for  beryllium  and  aluminum  was  about  100°C.  This 
comparatively  small  temperature  decrease  is  caused  by  the  high  thermal 
conductivity  of  metals  and  small  size  of  the  sphere.  As  a  result  of  the 
small  temperature  lag,  little  error  will  be  encountered  if  the  metal 
sphere  is  regarded  as  a  heat  sink  with  infinite  thermal  conductivity. 
Under  this  assumption,  the  rate  of  temperature  rise  of  the  surface  of 
the  metal  sphere  T  is  equal  to  the  total  heat  transfer  rate  to  the 
sphere  divided  by  the  total  heat  capacity: 


dT 
_ m 

dt 


6  k°  (T  1/2  +  T  1/2)  (T  -  T  )  + 

_ f  Jig _ m  g  °» 


oe  T  -  ae  T 
m  g _ m  m 


D  C 
0  P. 


m 


m 


(14) 


From  equation  (14)  a  computation  of  the  elapsed  time  to  heat  the  metal, 
T^,  from  Tq  to  T^  was  made  with  a  digital  computer.  In  doing  so,  the 
Dulong-Petit  limiting  heat  capacity  for  solids  was  used  to  compute  aver¬ 
age  specific  heats  of  the  metals  for  the  temperature  range  Tq  to  T^. 

The  additional  heat  capacity  of  liquid  metals  should  be  included  at 
temperatures  above  the  metal  melting  point.  However,  the  increase  is 
usually  only  about  1  cal/mole  or  157.  and  will  be  partially  offset  by  the 
lower  heat  capacity  between  the  initial  temperature  Tq  and  the  Debye 
temperature . 
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Delay  Due  to  Latent  Heat  of  Fusion 


An  additional  time  must  be  added  for  latent  heats  of  fusion 
of  the  metals  which  melt  at  or  below  their  respective  ignition  points. 


AHr  D  p 
fus  o  m 


(qv  +  qf>  t 


t  -  m:p.6 

m 


(15) 


The  computed  elapsed  times  for  metal  sphere  heating  and  fusion  are  also 
listed  in  Table  II.  The  sum  of  these  is  the  ignition  delay.  The 
temperature -time  relations  for  this  period  are  plotted  in  Figure  4. 


Effect  of  Particle  Diameter 

As  shown  in  Table  II,  the  heat  transfer  contributed  by  radi¬ 
ation  is  only  about  1%  or  27,  of  the  total.  If  radiation  is  neglected, 
the  elapsed  time  required  to  heat  the  particle  to  the  ignition  point 
under  stagnant  gas  conditions  is  proportional  to  the  square  of  the 
diameter.  The  heat  required  to  reach  T^  is: 


m 


(T, 


T  )  + 
o 


AH 


fus 


(16) 
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V. 


23 


FIGURE  4.  COMPUTED  IGNITION  DELAY  IN  A  METALLIZED  FLAME 


substituting  equation  (6)  in  equation  (16)  yields  the  parabolic  relation, 


t  » 


C  p  (T  -  T  )  +  AH. 
p  m  '  i  o  fus 

m _ 


12  k.  (T  -  T  ) 
f  g  m 


(17) 


where  (T^  -  T^)  is  an  effective  average  temperature  gradient. 


Comparison  with  Experimental  Results 

A  parabolic  dependence  of  the  time  on  ignition  particle  size 
has  been  observed  by  Friedman  and  Macek^  who  measured  the  ignition  delay 
for  aluminum  spheres  of  various  sizes  using  a  combustion  torch  similar 
to  a  Meeker  burner.  The  fuel  was  a  propane  and  oxygen  mixture  with  a 
calculated  adiabatic  flame  temperature  of  2510 °K.  For  a  100  micron 
sphere,  the  ignition  delay  was  21  milliseconds  (extrapolated). 

A  computation  of  the  ignition  delay  for  these  conditions  was 
made  using  equation  (14)..  The  thermal  conductivity  was  calculated  in 
the  same  manner  but  based  on  a  gas  consisting  of  57  mole  7«  H^O  and  43 
mole  7>  CC>2  which  reflects  combustion  of  propane.  The  lower  thermal 
conductivity,  caused  by  the  absence  of  low  molecular  weight  combined 
with  the  lower  flame  temperature,  2510°K,  extend  the  ignition  time.  The 
computed  result  was  8.3  milliseconds. 

Although  ignition  delays  were  not  measured  in  this  laboratory, 
the  maximum  aluminum  particle  size  that  ignited,  115  ±  10  microns,  was 


compared  with  the  maximum  exposure  time.  The  latter  was  calculated  from 
the  particle  velocity,  using  motion  pictures,  and  maximum  path  length  of 
the  particle  in  the  hot  zone  of  the  flame.  The  exposure  time  for  a  115 
micron  aluminum  particle  varied  between  15  and  20  milliseconds.  The 
computed  ignition  delay  for  this  size  was  11  milliseconds. 

When  the  hazards  of  estimating  thermal  conductivities  of  gases 
at  high  temperatures  are  considered  the  agreement  between  observed  and 
computed  ignition  delay  is  probably  satisfactory.  The  minimum  flame 
temperature  that  would  cause  ignition  of  aluminum  particles  was  deter¬ 
mined  by  Friedman  and  Macek  and  by  the  authors.  In  both  cases,  the 
minimum  flame  temperature  was  observed  to  be  within  ±  80°C  of  the 


melting  point  of  Al^Oy 


COMBUSTION  AT  ATMOSPHERIC  PRESSURE  WITHOUT  METAL  EVAPORATION 


It  was  previously  mentioned  that  aluminum  has  a  normal  boiling 
point  of  2740 °K  while  its  oxide  melts  at  a  lower  temperature,  2300°K. 

The  heat  flux  is  expected  to  increase  considerably  at  2300°K  because  the 
thermal  flux  contribution  from  the  exothermic  oxidation  reaction,  which 
is  negligible  below  2300 °K,  becomes  a  major  heat  source.  With  the  addi¬ 
tion  of  this  heat  source,  the  elapsed  time  analysis  is  similar  to  that 
for  the  ignition  delay. 

The  oxidation  rate  will  be  limited  by  diffusion  of  reactants, 
either  oxygen  or  aluminum,  through  the  oxide  layer.  Specific  information 
on  diffusion  in  liquid  A^O^  not  Available.  The  diffusion  coefficient 
varies  exponentially  with  temperature,  but  for  liquids  the  activation 
energy  and  argument  of  the  exponent  are  small  and  the  temperature  effect 
is  usually  slight.  As  a  result,  virtually  all  liquids  which  do  not 
possess  abnormally  high  viscosities  (glasses,  polymers)  have  remarkably 
similar  diffusion  coefficients  regardless  of  the  temperature  or  chemical 

composition.  The  value  of  the  diffusion  coefficients  are  usually  10 

_  in-5  2  -1  17 

to  10  cm  sec  , 
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The  oxidation  rate  and  associated  thermal  flux  of  aluminum 


decrease  with  increasing  thickness  of  the  liquid  oxide.  For  thin  oxide 

layers ,  (R  -  R  )  <  <R  , 
o  m  o 


dm /A 
dt 


R  )  — 
m  p  t 


(18) 


The  concentration  gradient,  AC,  depends  only  on  the  pressure,  which  is 
constant,  and  temperature.  The  oxidation  rate  may  also  be  considered 
in  relation  to  the  spherical  geometry  and  regression  rate  of  the  metal 
radius . 


,  /  .  un 

dm/ A  _ m  p 

dt  =  dt  m 


(19) 


Substituting  equation  (19)  in  equation  (18)  yields 


dr 

dt 


-d12ac 


p2  (R 
m  o 


which  upon  integrating  gives 


pi  d12ac  1 


(20) 


2 


2 


+  Const 


(21) 


and  installing  the  boundary  conditions,  at  t  °  0,  R  *  R  yields: 

mo 


(R 


R  ) 

m 


[ 


2  P  D  AC 


-J 


1/2 


(22) 


m 


Provided  (R  -  R  )<<R  ,  the  thermal  input  required  to  heat 
o  3 

the  metal  is  (  AH_  -AH-  )  -  «R  p  .  The  chemical  heat  source  is 

\l  \  3  o  " 

obtained  from  equation  (22)  and  consideration  of  the  particle  geometry: 


‘R 


“Ro  [2piB12 


AC  t 


1/2 


ah. 


(23) 


where  AH  is  the  average  heat  of  formation  over  the  temperature  range 

T.  +  T. 


BP 


(Tgp  -  T^)  .  For  ease  of  calculation,  the  value  of  (q^  +  q^)  at  2 
was  used  for  the  convective  and  radiant  heat  fluxes.  The  equation  for  the 
elapsed  time  for  this  period  is: 


4  ttR 


(AH  -  AH  ) 
BP  Ai 


—  p  ■  4  Jt  R^ 
m  o 


2  Pt  D12AC  t 


1/2 


AHf  + 


2 

4  R  t 
o 


(q  v +  qf} 


T  +  T 
i  BP 


(24) 


The  elapsed  times  calculated  from  equation  (24)  and  the  equivalent  values 

of  the  radius  regression  (R  -  R  )  are  listed  in  Table  III.  Calculations 

o  m 

were  made  for  three  values  of  D^AC;  10  10  ^  and  10  ^  g  cm  *  sec 

The  periods  are  somewhat  smaller  than  those  required  for  the  ignition 
delay  period  in  aluminum. 
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TABLE  III 


COMBUSTION  OF  ALUMINUM  BETWEEN  THE  MELTING  POINT 
OF  A1203  (2330°K)  AND  THE  METAL  BOILING  POINT  (2740°K) 


AC  (g  cm  *  rec  *) 


-4 

-5 

-6 

Property 

10 

10  J 

10 

AHf 

6400  cal/g 

6A00  cal/g 

6400  cal/g 

(^sp  * 

115  cal/g 

115  cal/g 

115  cal/g 

Elapsed  Time 

(milliseconds) 

.016 

.085 

.19 

(R  -  R  ) 
o  m 

(microns) 

.332 

.241 

.115 

METAL  EVAPORATION  AND  EXPANSION  OF  THE  OXIDE  SHELL 


Analysis 

At  the  boiling  point,  metal  Is  evaporated  from  the  metal  sphere 

and  transferred  to  the  oxide  shell  where  oxidation  occurs  as  shown  in 

Figure  5.  The  oxide  shell  must  expand  until  the  mass  diffusion  rate  of 

reactant  through  the  oxide  is  equal  to  the  equivalent  mass  evaporation 

rate.  If  the  evaporation  rate  decreases,  surface  tension  will  force  the 

oxide  shell  to  contract.  Metal  evaporates  from  the  surface  of  the  metal 

sphere  at  the  boiling  point  T  and  diffuses  to  the  outer  oxide  shell  at 

Ur 

a  higher  temperature  T£.  The  heat  removal  by  the  sum  of  evaporation  and 
transpiration  of  the  metal  vapor  must  maintain  the  sphere  at  Tfip.  The 
mass  transfer  rate  of  the  metal  can  be  expressed  in  terms  of  the 
radial  metal  vapor  diffusion  velocity 


W 

r 


dm 

dt 


4  ,t  r^  p  V 
v  r 


(25) 


\ 

\ 

»  >»• 
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FIGURE  5.  MODEL  FOR  COMBUSTION  WITH  METAL  EVAPORATION  AND  EXPANSION 
OF  THE  LIQUID  OXIDE  SHELL 
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It  is  assumed  that  is  small  enough  and  the  oxide  shell  sufficiently 
flexible  that  no  pressure  gradient  occurs  between  the  surfaces .  Then 
the  differential  equation  describing  the  temperature  distribution  T(r) 
in  the  space  between  metal  sphere  and  oxide  shell  is 


'dT 

av 

d 

1  r2  dT 

dr 

2 

r 

dr 

1  dr 

(26) 


The  thermal  conductivity  of  the  metal  vapor  is  assumed  constant  and 


calculated  at 


T  +  T 
BP  c 


from  the  theory  of  monatomic  dilute  gases. 


By  substitution  of  equation  (25)  in  (26) 
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W  CD 
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dr 


r2  dT 
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By  integrating  twice  and  using  the  known  boundary  temperatures 


we  find: 

I  -  Tc  __  e  ~kl/r  _  e  *kl/RX 

TBP  "  Tc  "kl/KR1  'kl/Rx 

e  -  e  ‘ 


(28) 


in  which  k^ 
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The  rate  of  heat  flowing  to  the  sphere,  Qv>  required  to  balance  the 
evaporation  heat  loss  is 


Q  =  W  AH  =  -4  it  K2IL 2  q_ 
v  r  vap  1  ^v 


r  =  KR, 


Since  for  spheres  the  heat  flux  is, 

.  dT 
%  v  dr 


then  by  substitution  in  equation  (29) : 


W  AH  -  4  jtk  K2  R  2  — 
r  vap  v  1  dr 


r  -  KR, 


From  equation  (28)  at  r  »  KR^, 


"kl/KR  "kl/R 

tbp  -  Tc  *  e 


From  this,  —  is  evaluated  at  r  »  KR^  and  substituted  in  equation  (30) 


to  yield 


Cp  <Tc-TBP> 
_ _v _ 

r fw  cn  (i  -  k) ■ 

J  r  Pv 

\  4  R  KR, 
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Equation  (32)  will  provide  in  terms  of  the  inner  oxide  radius  and  K 
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r 


(1  '  K) 


KR, 


4  jr  k 
_ v 
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<Tc  -  TBP> 


AH 

vap 


+  1 


(33) 


It  is  possible  to  determine  R2  f  (R^,  K)  from  the  mass  burned  and 
spherical  geometry 

Rq3  -  (KR^3  Pm  =  (R23  -  Rx2)  P]L  (34) 


For  small  values  of  AC,  the  equivalent  rate  of  diffusion  through  the 
spherical  oxide  shell  also  gives  the  mass  transfer  and  oxidation  rates 


4  «  D,  „  AC 


(35) 


The  unique  values  of  K,  R^,  R2  and  exist  where  the  mass 
transfer  rates,  from  equations  (33)  and  (35),  are  equal.  A  solution 
was  made  using  a  digital  computer  and  Newton's  method  of  solving  the 
following  equation  for  K  based  on  equations  (33)  to  (35) . 
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The  rate  of  heat  generation  by  reaction  at  the  oxide  must  equal 
the  rate  of  heat  transferred  to  the  surroundings  and  to  the  metal  sphere. 

Because  of  the  heat  capacity  of  the  oxide,  the  heat  of  reaction 
AH^  decreases  linearly  with  increasing  temperature.  This  effect  is  small 
and  calculations  using  reasonable  estimates  of  the  various  oxide  radii 
(Rq  »  50  microns)  would  require  an  oxide  temperature  T£  greater  than  the 
adiabatic  flame  temperature  for  aluminum,  magnesium  and  beryllium  if  only 
this  effect  is  considered.  This  means  that  AH^  must  be  limited  by  a 
shift  in  the  oxidation  equilibrium  or  by  partial  vaporization  of  the 
oxide.  The  equilibrium  constants  for  either  type  of  reaction  vary 
exponentially  with  temperature. 

Once  a  temperature  is  reached  where  these  effects  cease  to  be 
neglibible,  AH^  should  decrease  rapidly  in  one  or  more  steps  to  zero  at 
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the  adiabatic  flame  temperature  for  the  metal.  In  lieu  of  more  precise 
data,  the  oxide  temperature  is  expected  to  approach  the  adiabatic 
flame  temperature,  T^,  for  highly  exothermic  metal-oxygen  systems.  In 
the  following  calculations  for  each  metal  system,  an  arbitrary  selection 
of  the  oxide  temperature  was  made:  T  =  T.  -  200°C. 

C  A 

The  evaluation  is  repeated  for  all  selected  values  of  R^.  The 

resultant  oxidation  rate,  W  ,  and  associated  values  of  R  were  used  to 

*  r  m 

determine  the  regression  of  the  metal  sphere  radius  versus  time  using 
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p  4  ix  R 
m  m 


-dR 
_ m 

dt 


(37) 


Results 

The  results  were  machine  computed  and  are  plotted  in  Figures 

6  and  7  for  aluminum  and  beryllium  respectively.  The  three  curves  in 

-4  -5 

each  figure  represent  three  estimated  values  of  AC;  10  ,  10  and 

10  **  g  sec  1  cm  ^ .  The  computed  elapsed  times  for  combustion  of  100 
micron  spheres  are  compared  with  the  results  from  torch  experiments  in 
Table  IV.  The  agreement  with  the  estimated  values  of  D^AC  i8  excel¬ 
lent.  The  agreement  confirms  the  selection  of  these  orders  of  magnitude 
values  for  AC.  If  boiling  points  were  reliably  known  at  actual 
chamber  pressures,  an  identical  calculation  using  the  proper  boiling 

point  for  T  is  expected  to  provide  the  elapsed  time  for  combustion. 

BF 
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FIGURE  6.  REGRESSION  OF  RADIUS  DURING  COMBUSTION  OF  100-MICRON 
SPHERES  OF  ALUMINUM  USING  EVAPORATION  MODEL;  SYSTEM 
PRESSURE  ONE  ATMOSPHERE:  T.  »  3700°K 


39 


FIGURE  7.  REGRESSION  OF  RADIUS  DURING  COMBUSTION  OF  LOO-MICRON 
SPHERES  OF  BERYLLIUM  USING  EVAPORATION  MODEL;  SYSTEM 
PRESSURE  ONE  ATMOSPHERE;  TA  -  3900°K 


TABLE  IV 


V 


COMPARISON  OF  EXPERIMENTAL  AND  CALCULATED 
COMBUSTION  TIME  (milliseconds)  100  |i ALUMINUM  SPHERE 


Experimental 

Calculated 

D12  AC  =  10"4 

d12ac  -  10"5 

d12ac  «  io'6 

35-60 

56 

117 

160 
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DISCUSSION 


The  calculation  of  the  elapsed  times  for  the  ignition  delay 
and  combustion  periods  requires  specific  knowledge  of  material  constants 
at  high  temperatures  which  are  not  available  and  would  be  difficult  to 
measure.  The  errors  in  estimating  these  values  are  contained  in  the 
final  elapsed  time  calculation. 

The  unknown  constant  with  the  greatest  effect  on  ignition 
delay  is  the  thermal  conductivity  of  the  gas  in  the  convective  boundary 
film.  The  calculated  results  agree  with  experimental  results  within  a 
factor  of  two  and  this  discrepancy  is  within  the  probable  error  in  the 
thermal  conductivity. 

The  calculation  of  the  combustion  period  is  more  hazardous, 

but  errors  in  the  diffusion  coefficient-concentration  difference  product, 

D12  AC  exceed  all  others  in  significance.  Errors  in  this  product  could 

conceivably  amount  to  one  or  more  orders  of  magnitude.  The  results  for 

-4  -5  -1  -1 

aluminum  indicate  AC  is  between  10  and  10  g  cm  sec  which 


seems  reasonable,  but  an  Independent  determination  of  this  property  Is 
needed  to  verify  the  method  of  calculation.  A  considerable  effort 
towards  determining  the  diffusion  coefficients  of  0=  and  the  cations 
in  fused  oxides  of  aluminum  and  beryllium  is  needed.  The  desired  in¬ 
formation  is  useful  in  other  fields  including  the  reduction  of  metals, 
where  transport  through  a  slag  phase  is  involved,  and  single  phase 
ceramic  technology. 
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